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Si(lll)-7x7 Surface 

A. AbdurixitQ, A. Baratoff and E. Meyer 
Department of Physics and Astronomy, University of Basel, Klingelbergstr. 82, 4056 Basel, Switzerland 

(Accepted in Applied Surface Science) 

Molecular dynamics simulations have been performed to understand true atomic resolution, which 
has been observed on the Si(lll)-7x7 surface by dynamic force microscopy in ultra high vac- 
uum(UHV). Stable atomic-scale contrast is reproduced in simulations at constant mean height above 
a critical tip-sample separation when monitoring the interaction force between tip and sample. Miss- 
ing or additional adatoms can be recognized in such scans, although they are less well resolved than 
native adatoms. The resonance frequency shift, as well as arbitrary scans, e.g. at constant force 
can be computed from a series of force-distance characteristics. By means of dynamic simulations 
we show how energy losses induced by interaction with an oscillating tip can be monitored and that 
they occur even in the non-contact range. 

PACS: 07.79.Lh; 31.15.Qg; 34.20.Cf; 61.16.Ch; 68.35.Bs 
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I. INTRODUCTION 

Atomic force microscopy (AFM) has become an impor- 
tant method in fundamental and applied surface science 
able to obtain local information about surface structure 
and interactions, especially on nonconducting samples. 
In this method the surface is probed by measuring the 
interaction with a sharp tip or changes produced by that 
interaction. Atomic resolution can be obtained if the tip 
can approach the surface in the attractive regime with- 
out making contact. In such non-contact experiments it 
is difficult to measure this interaction force accurately 
while avoiding cantilever jump to contact. However, be- 
cause the resonance frequency of the cantilever changes 
locally if the tip starts to be attracted to specific sur- 
face sites, it became possible to image surface atoms by 
maintaining a constant frequency shift A/ via feedback 
control of the cantilever perpendicular displacement (!],§] . 
This shift can be measured very accurately using fre- 
quency demodulation Q . Jump to contact is avoided by 
using a large cantilever force constant k and a large tip 
oscillation amplitude A, such that the restoring force kA 
considerably exceeds the attractive force at the distance 
of closest approach. Under this condition, the tip-sample 
interaction weakly perturbs the harmonic motion of the 
tip, so that A/ can be calculated from the force-distance 
dependence at any particular site, as first pointed out by 
Giessibl @ 

Following the pioneering experiments many in- 

vestigators have obtained atomically-resolved images of 
the Si(lll)-7x7 surface in ultra high vacuum (UHV) 
using large-amplitude dynamic force microscopy (DFM) 
under various feedback conditions [l^] . In particular 
Liithi et al || recorded simultaneous images of the to- 
pography, of the time-averaged tunneling current and of 
the cantilever excitation, i.e. the damping, at constant 



A and A/. All images showed characteristic features 
(adatoms and corner holes) of the dimcr-adatom-stacking 
fault (DAS) model of the 7x7 reconstruction Point 
defects, e.g. missing adatoms appearing at the same 
locations in all three images; this convincingly proved 
true atomic resolution. Surprisingly, the contrast in the 
damping appeared inverted with respect to the other two 
images. Other issues are whether A/ jumps at the onset 
of covalent bonding between dangling bonds at the tip 
apex and on adatoms of the sample |fj"o|| , and whether 
optimum atomic-scale resolution is obtained at a closer 
separation |10|, in a soft intermittent contact mode 
or even in the range where A/ increases with decreasing 
tip sample separation ^3|. In this contribution we show 
how relevant information about those issues and quanti- 
ties measured in large-amplitude DFM can be effectively 
extracted from classical molecular dynamics simulations 
despite the enormous discrepancy in time scales between 
the tip oscillation and atomic motions. 

The rest of the paper is organized as follows: the next 
section describes the model and computational details; 
results are presented in the following two sections; con- 
clusions are summarized at the end. 



II. THE MODEL 

The Si(lll) sample is represented by a slab of eight 
layers with atoms coupled via the well-known short-range 
Stillinger- Weber (SW) potential Q which, besides a pair 
interaction, contains a three-body term favoring tetrahe- 
dral covalent bonding. Periodic boundary conditions are 
applied laterally to a nearly square supercell encompass- 
ing four 7x7 unit cells. Atoms in the two bottom layers 
are fixed in bulk-like positions. Atoms in the following 
three layers are allowed to move, but their average ki- 
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netic energy is controlled by a thermostat. The Si(lll)- 
7x7 reconstruction was initiated by removing and shift- 
ing atoms in the next two layers, and placing adatoms in 
the top layer so as to approximately obtain the geome- 
try of the dimer-adatom-stacking fault (DAS) model |p"l| . 
Dimcrs spontaneously formed between faulted and un- 
faulted half-cells in the third layer and adatoms equi- 
librated 2.7 'A away from their second layer neighbors 
shortly after molecular dynamics was started. Subse- 
quent simulations were initiated with this equilibrated 
structure. 

Focusing on short-range covalent interactions, which 
presumably give rise to atomic-scale contrast in the non- 
contact regime, we model the tip by a sharp pyramidal 
silicon cluster with 34 atoms in 6 layers. Like AFM tips 
etched out of Si, it has a [001] axis. They were kept at 
bulk-like positions in order to isolate effects due to sample 
deformation and dynamics. Because the SW potential 
vanishes beyond 3.8A , the longer-ranged pair part of a 
potential Jl3| which was fit to a first-principles computa- 
tion of the interaction between Si layers was added when 
considering the interaction between the tip and sample 
atoms. This potential has a decay length A=0.8A close 
to that found in similar recent calculations [[l4j of the in- 
teraction between a Siio cluster tip and a Si(lll) surface. 
Compared to Ref. 14, the force F z computed for the same 
Siio tip above an adatom using our interaction has a sim- 
ilar tail, but reaches a minimum (maximum attraction) 
at a distance z =3.0A from the unperturbed top layer 
which is too strong (-7 nN vs. -2.25 nN). Similar val- 
ues are obtained with our model tip. In the non-contact 
range, our results are thus at least semiquantitatively 
correct. 

The coordinates of the all moving atoms are up- 
dated according Newton's law, using a 5 t/l -order Gear 
predictor-corrector algorithm |15[| with a time step of 
0.4fs as advocated by Stillinger and Weber p"2fl . The 
average temperature is controlled by a first-order feed- 
back loop Jig ]. The thermalization time tt can be cho- 
sen such that the coupling to the thermostat is weak, but 
the discrete vibrational spectrum of our finite system is 
broadened without appreciable distortions. Because the 
SW potential led to appreciable surface diffusion at room 
temperature (which is not realistic), the setpoint temper- 
ature was fixed at 100K to ensure only small oscillations 
about the 7x7 equilibrium structure. Another impor- 
tant concern is the discrepancy between practical times 
for simulations and relevant ones in DFM experiments, 
e.g. the tip oscillation period 1// (few ps vs. few ms). 
Because A/ and the distance-dependent contribution to 
the damping Q arc induced by the tip-sample interac- 
tion, the cantilever dynamics per se can be ignored if 
only its stationary oscillation on resonance is considered. 

In order to track tip-induced changes, in each sim- 
ulation we recorded snapshots of atomic positions, as 
well as the time evolution of the following quantities: 
the center-of-mass of all moving atoms, their kinetic en- 
ergy (temperature), the coordinates of the tip apex and 



of selected adatoms, all components of forces on the tip 
and on those adatoms, the work done on the tip and the 
heat transfer to the thermostat. 



III. QUASISTATIC SCANS 



According to perturbation theory 
frequency shift independent of k can 
non-contact range from the z-dependence of F- 
ing possible reversible deformations), 



, a normalized 
oe computed in the 
(includ- 
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fo / F z [z(t)]cosujtf dt 
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where lj = 2nfo, z = d + A{1 + cosujt), d is the mini- 
mum tip-surface distance. In the large amplitude limit 
A > > d, constant A/ approximately corresponds to con- 
stant F z {d) at the turning point z — d. Quasistatic 
force-distance characteristics were obtained by approach- 
ing the tip at a constant velocity v z « v c — Xf ma x — 
lOOOm/s, f m ax = 16THz being the highest natural fre- 
quency of our system (adatom vibration against neigh- 
boring atoms). Such approach curves are essentially v z - 
independent and reversible down to a site-dependent sep- 
aration z c (e.g. 3.75A above adatoms) but strongly 
^-dependent and hysteretic below because irreversible 
symmetry-lowering jumps of the nearest adatom(s) are 
then induced, leading to a transverse force in the nN 
range. Although our interaction overemphasizes such in- 
stabilities, we expect them to develop below the inflection 
point in F z {z), even before the net attraction reaches a 
maximum. Indeed, this inflection signals the appearance 
of repulsive forces on the nearest adatoms which can be 
efficiently released by sideways displacement. 
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Long Diagonal [A] 



FIG. 1. Variation of the normal force on the tip along the 
long diagonal of the Si(lll)-7x7 unit cell recorded at 2 A/ps 
at a height z = 5 and 4 A measured from the unperturbed 
adatom layer. The symbols CH, A, R, X, D denote the posi- 
tions of the corner hole, corner and center adatoms, rest atom, 
hollow site and dimer row(short diagonal) of the underlying 
reconstructed surface 
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Below the highest z c constant height scans recorded 
at constant velocity exhibit friction and wear, adatoms 
being successively picked up and sometimes redeposited 
by the moving tip. Above that critical height such scans 
are reversible and, as shown in Fig. [I|, exhibit a con- 
trast with site-specific features in F z which becomes more 
pronounced at smaller z. Whereas maximum attraction 
above adatoms is expected fBLH, note that F z is the 
same above the rest atom (R) and above the hollow site 
(X) in the adjacent triangle between adatoms, although 
a recessed dangling bond exists only at the R-site. This 
implies that attraction to the adatoms dominates the 
picture. Because our model takes no account of charge 
transfer effects |Q , both cell halves and all adatoms ap- 
pear equivalent. 

Force-distance characteristics recorded at closely 
spaced points along a particular line can be used to re- 
construct scans measurable under all possible feedback 
operation modes. We checked that the constant height 
scan at z—A A in Fig. |] essentially coincides with the 
corresponding scan reconstructed in this fashion. The 
same procedure can be used to reconstruct scans at con- 
stant normalized frequency shift using Eq.(Q) with the 
time variable eliminated in favor of z = Acosujt in order 
to take advantage of the dense set of data points recorded 
in approach simulations, 
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The rhs of Eq.(||) is independent of spring constant k 
and oscillation frequency. From the reconstructed topo- 
graphical scans along the long diagonal shown in Fig. || 
we see that the constant frequency shift scan is close to 
the constant force curve. The contrast of the former is 
slightly reduced due to the non-uniform weight multi- 
plying F z in rhs of Eq. (Q) . As expected, the overall con- 
trast is inverted with respect to the constant height scans 
in Fig. [l| Furthermore, the apparent height difference 
between corner hole and adatom sites is approximately 
1.7A , i.e. larger than in measurements at constant A/, 
except one ]Kj . This is likely due to our neglect of inter- 
actions other than covalent bonding. The contrast we ob- 
serve is the maximum possible one for a given normalized 
frequency shift or force; the inclusion of Van der Waals 
and long range electrostatic interactions would reduce the 
contrast. The main advantage of constant height scans 
is that the contrast should be independent of such long- 
range interactions, because they are mainly determined 
by many atoms away from the tip apex. 

In order to assess the resolution of native surface de- 
fects, we first created representative ones by (i) removing 
a center adatom and (ii) placing an additional one at one 
of lowest energy surface sites H3 jlTj which is next to a 
corner adatom, and then relaxed the sample. Afterwards 
we recorded the constant height scans plotted in Fig. || 
together with the scan above the perfect surface already 
shown in Fig. [|. Compared to the native adatom sites, 



we observe an almost twice as strong attraction above 
the dimer formed between the additional adatom and the 
nearby corner adatom, but it is not possible to differen- 
tiate the two partners. On the other hand, the missing 
adatom produces a maximum almost as high and slightly 
broader than the corner hole. This appearance is consis- 
tent with that of deep minima associated with a missing 
adatom in constant A/ images M. 
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FIG. 2. Topographical scans at constant force (dashed 
curve) and constant normalized frequency shift (full curve) 
recorded under the same conditions as Fig. 1. Both 
curves were computed by interpolation from quasistatic 
force-distance characteristics. 
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FIG. 3. Constant height scans recorded under the same 
conditions as Fig. 1 in the presence of a missing center adatom 
(short-dashed curve, z=4 A), an additional adatom bound to 
a corner adatom (full curve, z=4.5 A) compared a similar scan 
above the perfect surface (dashed curve, z=4 A). The position 
of the missing (additional) adatom is indicated by the right 
(left) arrow on the horizontal axis. The oscillations appearing 
at the minima above adatoms are induced by the scan onset 
or by thermal fluctuations. 
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IV. PERIODIC TIP OSCILLATION 

In order to determine the distance-dependent damping 
attributed to the excitation of phonons || , we have per- 
formed dynamical simulations in which the tip was sinu- 
soidally driven so that it comes within a few angstroms 
of the sample. As expected, the time-dependent force 
on the tip consist of a periodic series of narrow pulses 
of width y/X/AJ f slightly modulated by thermal fluc- 
tuations. Fig. [| illustrates an extreme case where this 
minimum distance is smaller than z Cl so that the center 
adatom under the tip would be destabilized in a qua- 
sistatic scan. In view of the(unrealistically) narrow pulse 
width, a jump to an energetically more favored position 
on the side of the original site occurs only after the ninth 
pulse. Under typical experimental conditions the pulse 
width would be w 1/is i.e. encompassing more than one 
million adatom vibrations, a quasistatic situation. Nev- 
ertheless, a small amount of energy can be transferred 
to acoustic phonons, either directly or via anharmonic 
decay of local vibrations excited by the force pulses. In 
our simulations interaction-induced damping due to both 
processess || can be obtained if the inverse interaction 
time j \J Aj\ lies between f max and the frequency of the 
lowest mechanical resonance of our slab (about 1 THz). 
Thus that the imposed frequency cannot be much below 
that lower limit. 

To visualize energy transfer we calculate the work done 
by the tip as W = — J F z v z dt. From Fig. || we see 
that apart from quasi-reversible spikes accompanying the 
force pulses, W exhibits an irreversible increase, which 
comes about because sample atoms do not quite respond 
adiabatically. This results in a slight asymmetry of the 
force pulses. This irreversible part of the work goes into 
vibrational energy which is then shared between differ- 
ent modes i.e. turns into heat. This transformation is 
somewhat assisted by our weakly coupled feedback loop 
which regulates the total kinetic energy of all thermalized 
atoms according to the equipartition rule. 

The thermostat time constant tt can be chosen shorter 
than the period so that equilibration is practically 
achieved between successive pulses. The heat transferred 
to the thermostat then looks like a staircase as a func- 
tion of time, and the average transfer rate coincides with 
the mean power transferred by the tip (overall energy 
conservation) as demonstrated elsewhere |Q. In those 
earlier simulations, we observed that for a given force 
at the minimum distance d the energy transfer rate was 
stronger at the CH site than at the X site, in appar- 
ent agreement with the experimentally observed inverted 
contrast in the damping. As expected the damping con- 
trast was in accordance with the topography if the tip 
was driven with the same amplitude, but coming within 
the same minimum distance d=3.4 A. The observation of 
a contrast in the damping [|| at a constant force F z (d) is a 
delicate issue because such a contrast would vanish if the 
substrate were treated as a continuum subject to delta 



function like pulses |19| . In recent dynamic simulations 
we managed to detect a damping roughly proportional to 
F z (d) in the non-contact range, as expected for a linear 
response. We found that it is maximum above adatoms, 
i.e. at sites where the distribution of F z among nearby 
atoms is most localized laterally. 



-3 
-4 
-5 
-6 
-7 



m 



n 



n 



n 



(a) 



10 20 30 40 
simulation time [ps] 



50 




-100 



50 



10 20 30 40 
simulation time [ps] 

FIG. 4. Top: Time-dependent normal force on the tip os- 
cillated at 0.4 Thz with an amplitude A = 20 A above a 
center adatom to a minimum distance 3.15 A from the un- 
perturbed adatom layer. Bottom: Corresponding evolution 
of the work on the tip; the energy transferred to the sample 
after each pulse is converted in to heat(see text). 



V. CONCLUSIONS AND ACKNOWLEDGMENTS 

Molecular dynamics simulations can be conducted and 
analyzed effectively to provide valuable information un- 
available in experiments, as well as to compute measur- 
able properties (resonance frequency shift and damping) 
and corresponding scan lines. This has been illustrated 
for the Si(lll)-7x7 surface. Bearing in mind that our 
model tip is extremely sharp and that long-range in- 
teractions are ignored, the agreement with experimen- 
tal constant A/ scans is quite satisfactory. As with any 
simulations, some details are sensitive to the underlying 
model. The interactions we assumed tend to unduly fa- 
vor adatom jumps along the surface and to the tip below 
a separation z c which appears somewhat large. Neverthe- 
less, such a critical distance is expected to appear before 
the net force reaches a minimum (maximum attraction) 
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upon approach. It is therefore still difficult to under- 
stand how stable DFM images at constant A/ < can 
be obtained following a jump in A/ flt| or in the range 
where A/ increases with decreasing tip-sample separa- 
tion . The inverted contrast sometimes observed in 
the topography Q and in the damping || also remains 
unexplained. New attempts to reproduce such measure- 
ments, with careful attention to possible feedback arti- 
facts and to the relevant range of setpoints are therefore 
very desirable. More work is also needed to understand 
the detailed dissipation pathways and to extrapolate the 
computed damping for comparison with measured values. 
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this simulation and his molecular dynamics code, and M. 
Bammerlin, R. Bennewitz, M. Guggisberg and R. Liithi 
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